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IHTRODUCTIOH 


Radiance recorded by any reaote tenting instrument will contain noiae which 
will coneiet of both eveteaatic and randoa variations. Systeaatic variations 
aay be due to sun-target-sensor geoaetry (e.g. Duggin 1985a, Rirchner and 
Schnetzler 1981) ataospheric conditions (e.g. Dave 1978, 1980, Schnetzler 1981) 
and the interaction of the spectral characteristics of the sensor with those of 
upwelllng radiance (e.g. Slater 1979, Duggin, Slater and Soaers 1980, Duggin 
1985 a.b, Markham and Barker, 1985). Randoa variations in the data aay be 
caused by variations in the nature and in the heterogeneity of the ground cover 
(e.g. Daughtry, Vanderbilt and Pollara 1982, Duggin 1974, 1983, 1985 a,b), by 
variations in ataospheric transmission, and by the interaction of these 
variations with the sensing device (e.g. Duggin 1985a). 

It is important to be aware of the extent of random and systematic errors 
in recorded radiance data across ostensibly uniform ground areas in order to 
assess the impact on quantitative image analysis procedures for both the single 
date and the multidate cases. It has been shown that random variations in 
irradiance and in reflectance characteristics (caused. for example, by 
variations in the nature and in the heterogeneity of ground cover) can cause 
variations in the discriminabilitv of vegetation stress (Duggin 1983) and that 
random variations in unresolved (sub-pixel sized) cloud can affect 
discriminability of agricultural targets (Duggin, Schoch, Cunia and Piwinski 
1984). Duggin and Schoch (1984) and Wardley (1984) showed that the impact of 
random variations in irradiance, ground reflectance and atmospheric 
transmittance on target discriminability can be angle-dependent. Systematic 
variations in radiance due to scan angle have been observed by many workers in, 
for example, even MSS data with a scan angle range of 11.56 (Kaneko and Engvall 
1977) and in AVHRR data which has a much larger scan angular range of 55 (e.g. 
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Dustin and Saunders 1964. Dustin and Piwinskl 1964). Tha causa of fchs antular 
dapandancs of spectral radianca (and tharafora of discrimlnabillty) is tha 
systematic variation in tha raflactanca properties of fround covar with 
illumination and with viewing angles (e.g. Bauar, at al. 1979, Coulaon 1966, 
Coulson at al. 1965, Dustin 1977, Esbert and Ulaby 1972, Kollenkark at al. 1962, 
Smith 1979, 1963. Suits and Safir 1972). In tha casa of amittad radianca, thara 
is a dapandancs of amissivity on view anslas (a.g. Jackson 1961, Kimes, at al. 
1960. Kimss and Kirchnsr 1983, Kimas 1981a. b 1963). Atmospharic scattaring and 
transmission also vary with viawing and with illumination angles (a.g. Turner 
1978. Dave 1978). Tha combination of these systematic variations in factors 
controlling radianca levels gives rise to upwelling radiance which varies with 
viawing geometry in a target-dependent manner (e.g. Kirchner and Schnetzler 
1981, Schnezler 1981, Duggin, 1985, Duggin, Lindsay and Sakhavat 1985, Duggin, 
Sakhavat and Lindsay 1985). There is the possibility that systematic effects 
may be corrected for if they are properly understood. 

It has been reported that there is a sensor dependent variation of about 2X 
in radiance levels across the TK scan line in a manner which is scan direction- 
dependent (iialila. et al 1984, Metzler and Malila 1965, Kieffer, et al 1985, 
Murphy, et al 1965). This effect will compound the effect of the scan angle 
dependence of sensor output which is due to the goniometric anisotropy of the 
scene radiance. The angle-dependent upwelling scene radiance may interact with 
both the spectral and spatial characteristics of the sensor in an as-yet 
undetermined manner (e.g. Duggin 1985a, b): we do not yet know whether this is a 
first order or a higher order effect. 

It was our intention in this study to examine the systematic and the random 
variations in digital radiance data recorded in each band by the thematic mapper 
over crop areas which were ostensibly uniform and which were free from visible 
cloud. The thematic mappers on Landsats 4 and 5 have narrower bandpasses and a 
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wider rant* of waveband* than the HSS or the AVHRR and ao finding* for the MSS 
and for the AVHRR cannot automatically be assumed to apply to the TO, even 
though the ecan angle range for the AVHRR be restricted to that of the TO. For 
example, the superior spatial resolution of the TO (30 m as compared to 79 m for 
the HSS and 1 km for the AVHRR) will Increase the random variation between 
radiance values recorded from individual pixels located in apparently uniform 
areas (e.g. Daughtry, Vanderbilt and Pollara 1982). 

The analysis was performed on several scenes at different growth stages. 
We considered agricultural areas. We wished to see if there were seasonal 
effects upon both random and systematic variations in digital radiance data 
recorded in the thematic mapper bandpasses. 

The heterogeneity of the scene within the ground resolution cell, which is 
imaged by the IFOV will affect the sensor output in each bandpass and will 
affect the spectral distribution of sensor outputs for each pixel. This is 
because of the spatial distribution of different scene elements with different 
optical properties and anisotropies within the IFOV. The scene elements will 
interact with the rear-projected point spread function of the TO system onto the 
object plane in a manner which is both spatially and spectrally-dependent. This 
has been discussed by (Duggin 1985a, b) and will be dealt with in greater detail 
later in this report. The variation in scene type, and in the distribution of 
scene components across the area imaged by the TO will superimpose random 
variations in the level and in the spectral (be tween-band) radiance 
distribution. Such variations in scene composition across the imaged area will 
also give rise to edge effects. The impact of the edge effects will depend not 
only on the nature of the contrasting composition of adjoining ground resolution 
cells, but will also depend upon the angle through which the imaged area is 
viewed and upon the interaction of the upwelling radiance with the spectral and 
soatial characteristics of the TO sensors. 
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The edge effecte will impact the TH radiance data (calibrated digital 
counte) in a manner which la dependant upon the modulation tranafer function 
(HTF) of the TH in each bandpaaa for each aenaor (e.g. Schowengerdt , Archwametry 
and Vrlgley 1985. Cuahnie and Atkineon, 1985). Tha Fourier tranaforai of the HTF 
deecrlbee the inatantaneoua rear-projected point epread function (PSF) on the 
object plane, aa will be diacueeed later. 

The radiometric calibration of the TH data hae bean the focue of eeveral 
etudiee (e.g. Tilton. Harkham and Alford 1985, Deeachy, et al 1985, Schott and 
Volchok 1985. Hurphy, et al 1985, Singh 1985, Duggin 1985b). There haa aleo 
been eoneiderable effort to develop algorlthme to deetripe and to 
radiometrically correct the data after acquieition (e.g. Bernatein, et al 1964, 
Poroe and Peteraon 1985, Malaret, et al 1985, Wrigley, et al 1964, 1965). Some 
conai derat iona will be diacueeed here. 

There have been etudiee to evaluate the multivariate atatiatical analyaia 
of digital TH multiband, multitemporal data for improved land uae analyaia 
(e.g. Anuta, et al 1984, Forater, et al 1965). Studiea of the analyaia of 
eigenvectora of the principal componenta of multi temporal, multichannel TH data 
will be diacueeed in an appendix aa a potential meane for enhanced land uae 
diacrimination. 
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ANALYSIS 


In this analysis, we were constrained by data availability; tha data which 
wa hoped to obtain at various stajras ovar agricultural regions including ona in 
Iowa (path«27, row«31 on tha World Rafaranca System (WRS)) was not availabla at 
all of tha growth atagaa raquastad. Soma acquisitions wars cloudy and 
unforasaan circumstancas pravantad tha acquisition of othar acanas in tima to 
parform tha analysis for this report. Wa did analyze two imagas ovar a forastad 
araa, howavar. 

Tha imagas which we discuss in this raport ara rastrictad to three regions. 
In each case, we used radiometricallv corrected p-type CCT data (NASA 1983). 
WR5 path 27, row 31 covers a corn/soybean region in Iowa: in August 1982, only 4 
band data was available for this area. 

Firstly, in each case, the data were screened using tha Landsat Assessment 
System for cloud and for uniformity. That is, to ensure that tha scene did 
indeed consist entirely of crop areas for those regions examined. Only the 
roads between fields ware not vegetated for the agricultural areas studied. 

In the first analysis, a mask was generated. Three swaths were used across 
tha full image: each was 300 lines deep and started at lines 500. 1900 and 5000. 
Slices which were 16 pixels wide were taken in these swaths. The slices had 
starting pixel numbers 300 (bottom two swaths) 500, 1000, 2000, 3000, 4000, 
5000, 6000, 6500, 6700 (top two swaths). The offset is related to the Earth 
rotation correction of the Landsat image. In this manner, a non-biased analysis 
was performed by analyzing all of tha pixels in each slice over an apparently 
uniform, cloud-free agricultural region. Training within these regions on areas 
which appeared uniform on the image was not performed in this analysis, as it 
was considered that this would have resulted in bias deriving from 
unsubstantiated, a priori assumptions as to the nature of the target. Tha mean 
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digital counts for each tiles were calculatsd for aach bandpass. Tha digital 
counts Mara usad sines we wars intarastad in variations within 1 stages and 
eonsidarad that errors due to offset would not seriously affect our estistatas 
of coafflclants of variation (standard deviation divided by tha means C.V.) in 
radiance, except to altar tha C.V. by about 15* in tha case of TM band 6 (Barker 
1984, Barker at al 1964). At tha sane tine, tha variance and tha coefficient of 
variation ware calculated for each slice. Fig. 1 shows tha stean digital 
radiance values for tha first four bandpasses of tha August 02 Inage for path 
27, row 31. Only four bands ware available for analysis for this inaga. The 
naan radiance values are shown as a function of scan angle (starting pixel 
value for the slice) end of neen seen line for the swath fron which the slices , 
were taken. In this and in all subsequent figures, the northernnost swath 
(lower neen line number) is closest to the viewer, west is to the left (lowest 
pixel nunber) and eaut is to the right (highest pixel nunber). There was no 
attenpt to register inages in this investigation. 

It is seen that there is a significant systematic variation with scan angle 
before harvest, with a superimposed random variation of about 5*. The 
systematic variation is over 10 * between the edges and the center of the image 

and is apparently close to symmetric about nadir for TM band 1 of the pre- 
harvest image. The effect becomes more pronounced in band 2 and is almost 25* 

in band 3. There is a strongly assymmetric 30* change across the image for TM 

band 4 of the same image. There appears to be a general decrease in mean pixel 
radiance from south (mean scan line 5150) to north (mean scan line 650) in this 
image, coupled with some change in the apparent scan angle dependence. Also 
shown are the coefficients of variation (C.V.s) for the pixel radiance values 
contained within each slice of each single-band image. The C.V.s are also 
plotted as a function of pixel and mean scan line for the slices considered in 
the mask superimposed on the image. There is approximately a 50* variation in 
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C.V. about nadir for banda 1 through 3, with a superimposed random variation in 
tha C.V. and a systematic decreasing trand towards tha north portion of tha 
image (daeraaalnt moan aean lino). Band A (tha raflactod infrarad ration) shows 
sialnly random variation with a slitht monatonically lncroaslnt trand in C.V. 

t 

from wast to oast. It is notaworthy that tha C.V. is tanarally lass than 10* for 
band 1 and 15* for band 2 but rlsas to naarly 30* for band 3 (whoso digital 
valuas ara lowar than bands 1 or 2), falling back to lass than approximately 17* 
in band A. Thar a doas appear to be a general trand for tha C.V. to dacraaaa 
from south to north. 

Tha same region (path 27, row 31) was viewed again after harvest (October 
21. 1982). A color infrared rendition of this image on the interactive computer 
screen suggested that this area was mostly stubble. The mean digital counts for 
the seven band data are shown in Fig. 2 as a function of the same variables as 
for Fig. 1. They appear to show a general trend decreasing approximately 10* 
from west to east in the image. The reverse is the case in band 6, the thermal 
Infrared channel. There appears to be a slight decreasing trend in radiance 
values from south to north in the image and the noise (random variation) in 
digital radiance values appears to be approximately 10*. The digital values 
are lower after harvest, except in band 3. The coefficients of variation for 
this image are shown as a function of pixel and mean scan line in Fig. 3. They 
are all slightly higher than in the case of the pre-harvest image for the first 
four bands, are around 20* in band 5, 20* in band 7, but less than 5* in band 6. 
However, there is an apparently anomalous increase in C.V. at the far east side 
of the image, which was not readily explicable from image data of the slices 
examined on the interactive computer screen. The most obvious possibility would 
be a greater heterogeneity in ground cover at the eastern edge of the scene. 
For this image, an examination of the scene and of the analyzed slices in a 
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false color rendition on tho interactive computer ecreen suppested that patchee 
of viporous vepetation exiated in what appeared to be etubble or coil areae. 
The dletribution of theee ecene element* mipht have, for aome reaeon, been more 
heteropeneous towards the extreme eaet of the imape. 

The tame area (path 27, row 31) waa examined uainp a later (Landaat 5) 
acquisition obtained on Aupuet 15, 1984 (imape number 5016716293). For thle 
imape, in order to avoid allpht, localized cumulue cloud it waa neceeeery to 
atart the three awatha at line* 2072, 3900 and 4900. Mean dlpltal radiance 
value* are ahown for the teat area* (slice*) as a function of pixel and of mean 
scan line in Fip. 4. Bands 1-4 show weaker systematic trends than the Aupust 
02, 1982 imape of the same area. Random variation appears penerally to be of 
the order 5* in dipltal radiance values, while there is no obvious symmetry in 
the systematic component of variations in bands 1-3. Band 4 shows approximately 
e 20* decrease in radiance for the southern portion of the imape, but no such 
trend in the middle or for the northern repion. The dependence on mean scan 
line seems pronounced only for bands 5-7. Band 6 does show the same peneral 
increase to the east as for the October 1982 imape, analyzed in Fip. 2. 

Fip. 5 shows the coefficients of variation of the pixel radiance values for 
the windowed areas (slices) described by the overlay mask, plotted as a function 
of pixel and of mean line. The C.V.s are below 0.08 for band 1 and penerally 
below 0.10 for band 2, with a random variation of up to 30* and with only a 
slipht systematic decreasinp trend to the northeast. In band 3, the C.V. is 
0.25 at the west edpe of the imape, fallinp to 0.15 or less towards the east. 
The decrease is more pronounced in the north of the imape than is the case in 
the south. In band 4, the C.V. is penerally loss than 0.18, with a random 
variation of up to 20* and a slipht decrease from west to east in the south of 
the imape. The situation is similar in band 5. The thermal infrared band, TM6, 
shows C.V. values less than 0.04, which exhibit a peneral decrease in trend from 
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west to east In tho south of tho image, and a general dacraaaa fro* south to 
north, which is mors pronounced in ths oast of tho lmaga. Band 7 has higher 
C.V. valuos (up to 0.40). Thoro is a docroaso frost west to oast and a slight 
decreasing trend frost south to north. However, while the thersial IR band (6) 
shows the northeast region of the istage to exhibit the lowest variance, bands 4, 
5 and 7 (reflected-to mid-XR) show high variance. This may indicate a higher 
heterogeneity in growth stage in this region. 

The analysis on the August 1962 istage, when costpared to that performed on 
the August 1984 image suggests that the systematic variations across an image 
depend upon time. This may be related to the substantial non-uniform changes 
with Julian date in both the level and angular dependence of radiance recorded 
over the crop areas of the U.S. Great Plains by the AVHRR, as reported by (e.g.) 
Duggin and Piwinski (1964). Atmospheric changes and variations at ground level 
can occur between image acquisitions. 

Factors contributing to radiance changes across an image are atmospheric 
changes across t’ : imaged area, together with atmospheric scattering anisotropy 
and hemispherical-conical spectral reflectance anisotropy which is dependent on 
sun-target-sensor geometry, as mentioned earlier. However, while these effects 
will be substantial for a large scan angle range, covering a large area, such as 
the AVHRR (±55° ), one would expect these effects to be less for the TM, whose 
scan angle range is only i 7.7°. 

It has been noted (Duggin 1974, 1983) that the random variation to be 
expected in recorded radiance will arise partly from random variations in 
atmospheric transmission and partly from variations in irradiance: reported 
coefficients of variation are approximately 6X (Duggin 1974, 1983). It has also 
been reported (Duggin 1963) that ground reflectance measurements made at 80 m 
spacings in the MSS bandpasses show between 5X and 20* coefficient of variation. 
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Systematic recorded radiance dua to atmospheric scattering and bidirectional 
raflaetanea factort anisotropy might ba axoaetad to give risa to substantial 
scan antis dependence for larts scan antis ranges (such as + 5^ for tha AVHRR) . 
However. random variations might ba axoaetad to pradominata ovar systsmatlc 
variations for tha smallar scan angular rants (±7.7°) for tha thamatic mappar. 
It is also intarasting to nota that tha ranga of random variation bafora and 
•ft ar harvast in 1982 is not markedly different and that tha same range of 
variation aooaars to aoolv to tha August 1984 images. 

Ua now triad a different tachniaua. involving ’'training'’ within tha grid 
call areas described above. Tha same mask (ovarlav) was employed as in tha 
method described above in order to define slices 16 pixels wide by 300 lines 
daac. However, in this analysis, we used tha cursor whose size, shape and 

position could ba defined interactively to "train*' on areas considered to 
consist of tha most vigorous vegetation because of their red hue and saturation 
whan tha Interactive disolav was used with appropriate look up tables (LUT’s) in 
the following modes: 

(i) TH 2 • blue 

TH 3 « green bast for agricultural scene 

TH 4 * rad 
(ii) TH 4 * red 

bast for forest scene - 

TH 5 ■ green separates forest from grassland 

and agricultural areas 

TH 7 » blue 

Each "slice" on each "swath" was taken to ba a class, and training was performed 
on three vigorous vegetation cites in each class, or slice. In this analysis, 

we did make the assumption that what appeared to be "red" on the monitor was 
visorous vegetation, so that this was a biased analysis. 

The cursor mode war used on the LAS interactive analysis terminal (IAT) in 
order to enclose those areas identified (g priori ) bv means of their appearance 
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In • f«lM color mods. Tho control point of the examined area was obtained 
(both acreen coordinate* and coordinate* related to the master imago) from th* 
1AT control monitor, using the P0LYSIT2 program. While th* number of pixel* 
enclosed in each slice was always constant, this quantity varied in the cae* of 
individually located and positioned cursors. 

Th* analysis described above was used on the pre-harvest Fort Dodpe imapes 
(path 27, row 31, imape # 4001716261, 02 Aupust 1982). Three trainlnp sites on 
apparently viporous vegetation were selected in each slice and th* data for th* 
pooled sites in each sliea are presented her*. Fip. 6 shows th* variation with 
pixel number (scan anpls) of th* mean dipital radiance values for each training 
sits, in each bandpass. Swath X has a line start (first line of swath) of 500, 
swath Y of 2000 and swath Z of 5000. (loan pixel radiance values for each of tha 
two training areas within aach grid call, as well as th* mean of the means of 
th* two training areas or cells (for which th* best fit curve is shown) ar* 
displayed. Tha scan angle-dependence is larger than in th* case for tha 
previous analysis. Scan angle effects of up to 2 OX are seen in TH bands 1 to 3, 
but th* affects in TH band 4 are anisotropic about nadir, and exceed 20X. Th* 
reason for tha difference in the results of the two analyses is probably dua to 
the fact that th* method involving no training on apparently vegetated areas 
Includes areas which are lass vegetated, or areas whose leaf area index is lass 
than that in areas (subjectively) selected in this (interactive) analysis. When 
trainlnp is employed, an assumption is mad* that vegetation vigor is related to 
hue and to saturation perceived on th* screen of th* image analysis terminal. 
Training within the (18 x 300 pixel) grid call areas was also used for the seven 
band post-harvest data for th* same seen* (path 27, row 31 image ♦ 4009716273, 
21 October 1982). Th* data for th* X swath (line start 500) is shown only in 
Fig. 7. It is seen that there is far less variation than was tha case before 
harvest. This is to be expected, since th* scene consists mainly of stubble. 
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An examination of a poet harvest Indianapolis , Indiana scene (image • 
4010315505, 27 November , 1962) shoved similar results to the post harvest 
tfebster-Fort Dodge image (path 27, row 31 image # 4009716273). The plot of the 
swan digital counts for each of the grid cells as a function of start pixel is 
shown in Fig. 6, and exhibits the sane general behavior as shown for image # 
4009716273, the post-harvest image for path 27, row 31. The plot of the coeffi- 
cients of variation in digital counts for each bandpass, for each 16 pixel by 
300 line grid cell, displayed as a function of start pixel, is shown in Fig. 9. 

He examined an image of forest at Jamestown. We did not train within the 
grid cells. The image number was 4004315244, August 28, 1962. One swath was 
taken across the image at (approximately) line 5000. The swath was chosen so 
that the slices contained forest and minimal cloud. When the ’’slices” (16 x 300 
pixels each) were displayed in composite form on the IAT, a very strong angle- 
dependence was obvious to the eye. The coefficient of variation for the mean 
digital count in each band and for each slice was generally less than 15*. Fig. 
10 shows the large anisotropic scan angle-dependence of the mean digital counts: 
the effect is approximately 30* in the case of TH band 3. TH band 4 shows 
scatter. Fig. 11 shows the coefficient of variation (C.V.) for the digital 
radiance values, plotted as a function of scan angles. The scan angle- 
dependence of the C.V. is approximately a factor of two for several bands. The 
procedure of training on areas of apparently high vegetation was tried for the 
same area, imaged September 13. 1982 (path 17, row 31, image # 4005915251). 
Significant, but somewhat less scan angle dependence was found in this case, as 
is shown in Fig. 12 for the mean values. TH 4 again showed scatter, while TH 6 
showed anisotropic behavior. The points for the individual training areas were 
closer to the best-fit line for the pooled data than was the case for the post- 
harvest agricultural data. The scan angle-dependence of the coefficient of 
variation is shown in Fig. 13. 
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COMMENTS ON THE INTERCALIBRATION OF MULTI SENSOR. MULTI TEMP ORAL . MULTICHANNEL 

DIGITAL RADIANCE DATA 


In order to compare recorded radiance data obtained on the same date using 
different sensors, or on different dates using either the same sensor or a 
combination of sensors such as the Landsat A and 5 sensors TM A and TM 5. it is 
necessary to refer such data to a common datum. This point is discussed in 
Duggin 1965b. 


While calibration is readily made for a sensor in which the gain and offset 
have been adjusted during operation (e.g. U.S. Geological Survey 1979, 1984, 
Duggin 1961), the intercomoarison between sensors with different spectral 
responses (bandoassee) is more complex (e.g. Duggin 1980, 1981, 1985; Slater 
1979). This is because of the interaction of the spectral response of the 
sensor with the spectral radiance incident upon the sensor and because different 
sensors may record information in different parts of the spectrum, where the 
upwelling radiance has different values. This point is shown for a 
hvoothetetical situation in Fig 14. Here bandpasses A and B from instrument 1 
are taken to be approximately equal to bandpass C of instrument 2. The bandpass 
of a sensor is normally taken to be the wavelength region between the boundaries 
where the sensor has a 50* response. It is seen that the spectral regions 
between the half-power wavelength limits of bandpass A of instrument 1( A-^ , Ag^) 
and between the half-power wavelength limits of bandpass B of instrument 
H A<jg , Agg ) are not (in sum) the same as that encompassed by the half-power 
wavelength limits of bandoass C of instrument 2( A lc , A 2c ). The spectral band 
shapes of the sensor response curves also differ and so the calibrated, summed 
outputs of the two bandpasses (A and B) of instrument 1 are 
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while the calibrated output of instrument 2 is 
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Here a^. eg and a^ are the offsets of the sensors in bandpasses A, B and C for 
the two instruments. Also. 1/b^.l/bg end 1/b^are the gains of the sensors with 
bandpasses A, B and C. X^,q . Xg . q and X^.q are the uoper 2 ero power wavelengths 
for the sensors, while X^q. XgQ and X^q are the lower zero power wavelengths, 
respectively. 


If L(X ) changes, which could be as a result of different viewing and 
illumination geometry (e.g. Kollenkark, et al 1962. Duggin 1985a, Smith 1983) 
then the interaction oetween L(x) and HX) will vary in an instrument-dependent 
manner. Thus, intercalibrations between the outputs from bands A and B of 
instrument 1 and band C of instrument 2 depend upon the spectral scene radiance 
and may therefore change with view angle across even a homogeneous scene. The 
intercalibrations will also be affected by unresolved cloud, which will modify 
L(x). However, such intercomparisons are wavelength-specific. The calibrated 
outputs from channels A and B of instrument 1 will vary differently from each 
other and from that of channel C of instrument 2. Therefore, the intercompari- 
son must be target-dependent, view angle-dependent and atmosphere-dependent, as 
well as illumination angle dependent (in the ootical-ref lective region). The 
same principles of the above argument would apply if radiance data from two 
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basically similar sensors (each possessing ona similar band paai) on two dif- 
farant instruments wara balnf compared. 

Furthar, at pointad out by Forshaw, at al (1983), by Duggin (1985a) and by 
Duggin and Schoch (1983), tha output of a given aanaor dapandt not only upon tha 
natura of tha target and ita heterogeneity, but upon tha ralativa location of 
acana elements of diffarant optical propartiaa and aniaotropiaa in tha nominal 
inatantanaoua field of view (1F0V) of tha aanaor. Thia ie becauae of tha 
lntaraction of the point-spread function (PSF) of the aanaor with tha radiance 
recorded froai each of the different ecene elements in the IFOV. Thia situation 
is shown diagrammatically in Fig 15. The scene is considered to consist of three 
components A, B and C which give rise to different radiance levels at the sensor 
(bandpass r) whose point-spread function is shown. The PSF fall-off to zero is 
generally Gaussian and so radiance may be recorded from very bright scene 
elements beyond the nominal IFOV of the sensor. The sensor output for a pixel 
will therefore depend upon the position of the PSF peak with respect to the 
assemblage of scene elements. 

As pointed out by Duggin and Schoch (1983) for a cloud-free image 
consisting of n different reflecting elements, each occupying a fractional area 
a n of a pixel (IFOV), the recorded radiance in bandpass r is, for the optical- 
reflective region, given by 

r™(rr>^rr If i(*,y) 

+ fj<*,y) + . . . + g„<xy) A„(6,*,#',4',\)P n (xj)] r(f>',A)dydiJ d6d# + ',>',,A)lriA 

Lr{8\4’) * — — — — - 1 — - l - I ( 3 ) 

v 

/(X).dA 

*40 


15 


where 1(A) 1* the spectral response of the sensor and A^q.A^qi are the lower and 
upper zero power bandpaet limits of the eensor (e.g., bandpass A), respectively 

E(6t^«X) is the scalar global spectral radiance Incident on the pixel 
T(0'’t '.A> is the atmospheric transmission along the path from the pixel to the 

sensor 

R (e.4 ; 6 \ A ) is the spectral hemispherical-conical reflectance factor for 

n 

scene element n. which occupies fractional area a n of the 
heterogeneous pixel 

* 

P^fx.y) is the point-spread function at position x,y in the pixel 
g n (x,y) is a delta function which equals 1 if scene element n is present at 
(s,y), but is otherwise zero, so that for an IFOV of dimensions X,Y. 

rr g n (xj)dydx m o„; 


L (e,.*„ » 6 d>'..A ) i* the oath radiance for the sun-scattering center-sensor 

p 1 Z Z 1 

geometry defined by (6-j • ^ • 6 '<| • 4 '<| ) which will differ from the sun- 
target-sensor geometry < e • 4> » 6 ' * d> ' >• 

Clearly, the effect of scene element heterogeneity and disposition (with 
respect to PSF peak value for a given IFOV) on pixel radiance values and on 
between-band radiance distribution depends on illumination and on viewing 
geometry. For multi temporal analysis, this consideration is important in 
dete» lining the acceptable (fractional pixel) superposition accuracy of multi- 
date images, so that the above effect does not adversely affect the accuracy of 
classification procedures. 

In the case of unresolved (sub-pixel sized) cloud, equation (A) will be 
modified by the addition of a term 
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J^/CX) | J|" £* J‘ J * ^»<* ^> + «*<*.?> W#T*M PA**)+-- 

+ inky) R m {t##A'Jk) Pn(xy) + ... Lr(f'A') • + lA*#) U.(X^teo')! f(g>',XWydj tj d»d* + L,(fli>i;fl,*'i,X)j4fX 

wr/i 



where g^ deecribee the presence or absence of cloud at a location in the pixel 
and where R^fX) is the (nearly Lambertian) reflectance factor of cloud. 

Similar equations can be developed for the thermal infrared region of the 
spectrum, where scene radiance is determined by scene element spectral 
emissivity e(6'.4>'.X) instead of by the product of global spectral irradiance 
E( 6 , 4>; 6 ' ) and the hemispherical-conical spectral reflectance 
factor R(6 ; 6 ' ' ,X ). 

In summary, more calibration data is required for remote sensing instru- 
ments. This should consist of gain, offset, spectral response and point-spread 
function for each bandpass. This data should be available even where radiance 
data have been radiometrically corrected, so that it is possible to understand 
what has been done to the raw radiance data. This information will permit 
theoretical studies which will lead to the development of sensor intercalibra- 
tion procedures and to an understanding of the limitations inherent in such 
intercalibrations. Such studies will improve understanding of the effects of 
unresolved cloud, viewing and illumination geometry, scene composition and 
heterogeneity. Perhaps correction procedures for systematic effects (e.g. scan 
angle differences) may be developed. Procedures for the determination of mini- 
mum acceptable image superposition accuracy for the multi-date analysis of 
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iMftM with various ssnsor-scana combinations, with assasiblagas of scans 
•lsmants with diffarant optical prooartias and diffarsnt hatsroiraneltias will 
also naad to ba davalopad. 
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CflECUffilfittS 


This study cannot bs considsrsd sxhaustivs: indssd it is still in progress 
as tha aultidata data continuas to arrive. However, savaral conclusions ara 
suggested by this work. Vhlla soaa systematic trends in radiance values with 
scan anal* ware observed prior to harvest over a crop area in 1962, the same 
pattern was not repeated two years later. It appeared that the random variation 
in mean digital values recorded from 4800 pixel sample areas at regular 
intervals across an image in three swaths generally exceeded the systematic 
variations for the three images studied, and that the coefficients of variation 
were within those which might be expected to occur from prior measurements. The 
coefficients of variation of the digital values from the 4800 pixel areas 
selected as regular intervals across an image showed some scan angle dependence, 
but were more dependent upon bandpass than upon season or upon scan angle. The 
systematic effects did appear to be significant for a forest area. 

Random variations may affect image classification accuracy. Further, 
uncorrected systematic variations across and between images may impose 
restrictions on the level of classification accuracy which may reasonably be 
expected from automated classification of single date or multidate, multichannel 
digital thematic mapper data for the quantification and identification of 
terrestrial features in a non-photointerpretive fashion. It is therefore 
important to understand the restrictions which such variations inherent in the 
digital radiance data may place upon analyses. To this end, further work is 
needed in which further empirical studies of digital radiance data are used to 
determine optimum regimes of data acquisition and analyses for selected feature 
identification and quantification. 

Training on areas of vegetation, selected on the basis of their false-color 
rendition has been shown to produce Indications of random and systematic 
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variations in digital radianca data for aach band pass, within and batwaan 
images. Howavar, this method la lata rallabla than a atratlfiad sampling 
technique, alnea tha method of training biaaoa tha data, dua to assumed 
ralationshlpa batwaan tha iaiaga eharactariatlca and ground covar. which nay or 
nay not ba valid. 

Tha affacta of lntarcallbration batwaan tha TH aanaora on Landsat 4 and 
Landsat 5 hava baan discuaaad. Thasa ara problems which hava attractad tha 
attantion of othar workara aa wall (a.g. Murphy, at al 1965, Singh 1985, 
Malarat, at al 19"5, Slatar 1979, Markhr.a and Barkar, 1985). Howavar, thara ia 
a naad to ba awara of tha intaraction of tha spaetral rasponsa charactariatica 
of tha aanaor and tha spaetral charactariatica of tha upwalllng scana radianca 
(a.g. Duggin 1985a, b), which can give riae to (a.g.) target-dependant, sun 
angle-dependent striping. There is also a need to reduce data to a common datum 
level for meaningful intercomparison (as also remarked by Metzler and Malila 
1985). 

Tha use of principal components analysis for multitemporal TM data analysis 
is considered in an investigation performed partly under the auspices of this 
contract and described in the paper attached in appendix A. There is promise 
that the evaluation of eigenvectors for certain ground features may enhance 
discriminability. 

Other papers arising from this work are included in appendix 6. 
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FIGURE CAPTIONS 


Fif. 1 


Fir. 2 


Fir. 3 


Fir. 4 


Fir- 5 


Fir. 6 


The means and coaffielanta of variation of digital radianca data for 
ractanrular 4800 pixel sample araaa (allcoa) daacrlbad in tba text, 
plottad aa a function of acan anrla (pixal) and of mean acan llna for 
tha pra-harvaat imara *4001716261. path 27, row 31, 02 Auruat 1982. 
Only four banda of data wara aval labia. 

Tha aiaan dlrltal radianca valuaa for tha ractanrular 4800 pixel sample 
allcaa of tha naak described in tha taxt, plottad aa a function of 
acan anrla (pixel) and of naan scan line for tha post-harvest imara 
#4009716273, path 27, row 31. October 21, 1982. 

The coefficients of variation of dirital radiance data for the 
ractanrular slices plottad as a function of scan anrla and of mean 
scan line for the imara shown in Fir. 2. 

Tha naan dlrltal radianca values for the ractanrular sample slices of 
tha nask described in tha text, plotted as a function of scan anrla 
(pixel) and of mean scan line for the pre-harvest, Aurust 15, 1984 
Landsat 5 thematic maooer imara #5016716293, path 27, row 31. 

Tha coefficients of variation of dirital radianca data for ractanrular 
slices defined by the mask plottad as a function of scan anrla and of 
mean scan line for the imara shown in Fir* 4. 

Tha mean dirital radiance values for tralninr areas selected within 
ractanrular imara sample areas, definad by the mask described in the 
taxt. Hare the imara is that obtained by Landsat 4 on Aurust 02, 1982 
for path 27, row 31. "Rad 1" and "rad 2" refer to mean dirital 
radianca values for individual tralninr areas and "rad" refers to tha 
avarara of tha mean values for those tralninr areas within a 'ch imara 
'slice'. 'X' *Y' and *Z' refer to 300 line swaths taken across the 
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Fig. 


Fig. 


Fig. 


Fig. 


Fig. 

Fig. 


nortb, center and south of tho image. Tha Man digital radianea 

valuta for tha first four TH bands art plotted as a function of pixel 

for sach swath. 

7. Tha Man digital radianea valuas for apparantly vagstatsd training 

arsas ("rad 1" and "rad 2") and for avaragaa of tha Mans from thoaa 
training araaa ("rad") takan from within rectangular image sample 
araaa dafinad by tha mask dascribad in tha taxt. Tha hataroganaity 
within aach si lea is contl darabla aa tha image was obtaintd Octobar 
21, 1982 ovar a crop araa path 27, row 31, whara much of tha scans 

conslstad of stubbla. Saa Fig. 2 for comparison. 

8. Tha Man digital radianea valuas for tha ractangular slicas of tha 

mask dascribad in tha taxt, plottsd as a function of p.v*l and moan 

scan lint for tha post-harvest Novsmbar 27, 1982 Landsat 4 imago 

#4010315505, path 21, row 32. 

9. Tha coafficiants of variation of digital radianea data for 4800 pixsl 
ractangular sampla slicas dafinad by tha mask dascribad in tha taxt, 
plottad as a function of pixsl and of maan scan lina for tha image 
shown in Fig. 8. 

10. Tha maan digital radiance valuas for tha ractangular 4800 pixel slicas 
of tha 300 lina swath starting at lina 500 dascribad in tha text, 
plottad as a function of pixel and maan scan lina for tha August 28, 
1982 image of tha Jamestown forest araa (image #4004315244, path 17, 
row 31 on tha VRS system). 

11. As Fig. 10, but for tha 300 lina swath starting at line 2000. 

12. Tha maan digital radianea valuas for apparently vegetated training 

areas ("rad 1" and "rad 2"), takan from within ractangular image 
slicas. These ware taken across one swath (west to east) of an image, 
and plottad as a function of pixel. Tha image (#4005915251) was 


23 


obtained by Landaat 4 September 13, 1982 over a forested area, path 
17. row 31. 

Coefflclonta of variation of digital radiance value* in each bandpaea 
for the image area* deacribed in Pig. 12. 

Spectral reeponeea of two different remote earning device*, on* of 
which ha* two adjacent channel* similar in combined wavelength 
coverage to a single channel on the other device, ahown with the 
upwelllng spectral scene radiance. 

Interaction of the roar-projected point spread function of a sensor 
with a heterogeneous scene composed of elements whose optical 
properties (and therefore, scene radiance values) differ. 
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Previous work has suggested that seasonally varying reflectance properties 
are predictably related to radiance recorded by siultichannel remote sensing 
devices. Two images were obtained from thematic mappers on Landsats 4 and 5 
over the Washington, D.C. area during November 1982 and March 1984. These were 
registered and selected training areas containing different types of urban land 
use were examined, one area consisting sntirely of forest. Mean digital 
radiance values for each bandpass, in each image were examined and variances, 
standard deviations and covariances between bandpasses were calculated. We 
found that two bandpasses caused forested areas to stand out from other land use 
types, especially for the November 1982 image. In order to evaluate quanti- 
tatively the possible utility of principal components analysis in selected 
feature extraction, the eigenvectors were evaluated for principal axes rotations 
which rendered each selected land UBe type most separable from all other land 
use types. The evaluated eigenvectors were plotted as a function of land use 
type, whose order was decided by considering anticipated shadow component and by 
examination of the relative loadings indicative of vegetation for each of the 
principal components for the different features considered. The analysis was 
performed for each seven-band image separately and for the two combined images. 
We found that by combining the two images, we obtained more dramatic land use 
type seoaration. Conclusions have been drawn from this preliminary work sug- 
gesting directions for further study. Both British and U.S. image analysis 


systems were used. 


IHTRODUCTIQM 


Thematic mapper (TH) data have a nominal spatial resolution of 29m in all 
bandpasses except the thermal infrared band. The TH potentially has utility for 
measurements of urban land use and of vegetation confiscation within land use 
types (e.g. Quattrochl 1963, Haack 1963, Forster 1963, Jensen, et al 1963). 
Registered, multidate i maces of TH data each contain 7 bands of data. There is 
promise of superior feature identification and of urban forest identification 
using principal components analysis of TH digital multiband radiance data (e.g. , 
Bernstein, et al 1984) or other transformations (e.g., Crist and Cicone, 1984) 
to reduce the dimensionality of the data without reducing its Information con- 
tent. This approach has been discussed frequently (e.g., Stiteler 1979) and 
will facilitate data reduction. Such considerations are important in economical 
data analysis and in the representation of mapped features using images formed 
from only three primary colors. 

In this report, we discuss the results of an experiment to show that 
principal components analysis may be used to effectively distinguish between 
different urban land cover classes and that this may be done more effectively by 
using a combination o f images obtained on two different dates than by using 
either of the images individually. In order to quantify the superiority of 
using combined (super! nposed) images, we have chosen to evaluate the eigen- 
vectors for selected test areas, showing that the magnitudes of the eigenvectors 
are most different for the combination of images obtained on two dates then for 
either of the individual images, in the cases of those principal components 
explaining most of the variance in the data. 


ANALYSIS AMD DISCUSSION OF RESULTS 

In this preliminary experiment we utilized two registered TN images 
obtained over the Waehinirton, D.C. area (November 2nd, 1982 and March 24, 1984). 
The former image wae obtained by Landeat 4 (image no. 4010915140) and the latter 
wee obtained by Landeat 5 (image no. 5002315112). We regietered theae image* on 
the Landeat Aeeeeement Syetem at NASA, Goddard Space Flight Center. The euper- 
impoeed 512 pixel by 512 Line extract* were then output to tap* and eubeequently 
were interactively analyzed by using the GEMS interactive image analysis eytem 
located at the Royal Aircraft Establishment, Farnborough, England, and also by 
using the GEMS interactive image analysis system at ERSAC, Ltd. in Edinburgh, 
Scotland. The individual images are shown in Fig. 1. A linear stretch has been 
applied to each image. 

Four training sites were selected as "typical” and are shown outlined in 
the lower center image in Fig. 2. These consisted of "Forest", "Dense urban", 
"Downtown" and "Airport". The forest area has a 100* canopy cover, and is to 
the northwest of the city. The downtown area includes the park areas which 
contain few trees but extensive lawns, including those between the Lincoln 
Memorial and the House of Congress. The dense urban area has virtually no trees 
and is to the east of the downtown area and Washington National Airport is 
included as the fourth area, being a mix of concrete, buildings and grass, the 
last being less well maintained than that in the area between the Lincoln 
Memorial and the House of Congress. These four areas are annotated on the image 
at the bottom center of Fig. 2. The training sites contained the following 
numbers of pixels. 


Forest 

7480 

Dense urban 

5852 

Airport 

2448 

Downtown 

9246 


These sample sizes were sufficiently large to perform a reliable analysis. 


II* have shown ths Min radisne* valuas recorded for aaeh of tha four land 
usa categories usad in tha study In Tabla 1. In this table, we hava ordarad tha 
land usa cataitorias in tarns of dacraasint vegetation contaNt. TM band 4 (0.76- 
0.90 n) shows a dacraasint valus with dacraaaint vatatation for both lnatas, 
although tha downtown araa doas appaar to show approxinataly tha a ana band 4 
radianca valua as forast in March 1984. This would saan raasonabla, sines 
daciduous trass are not in laaf in March, and will tlva risa to lowar band and 
radianca valuas. TM band 6 (10.4-12.5 ») shows tha thsrnal radianca racordad 

fron tha training areas increasing with dacraasint vatatation contant. Again 
tha distinctions ara nost consistent in Novambar 1982. A ratio of (TM band 6/TM 
band 4) shows an incraasint value with dacraasint vatatation contant. This 
suggests that TM imagery could ba useful in mapping vatatation. especially in 
Slapping forast in urban environments. 

Tha principal components transformation was performed by using training 
areas (Fig. 2) over selected land usa types, so as to rotate tha measurement 
axes (i.e., axes describing recorded radianca in each band) to make tha selected 
land usa type most separable from "everything else". Tha eigenvectors ware 
evaluated for each of the four selected land use types identified in tha train- 
ing sites after principal components analysis was performed to emphasize that 
feature. Similar calculations wars performed, but without training, for the 
whole 512 pixel x 512 line image extract: the 14 principal component images for 

the two superimposed images are shown in Fig. 3. The eigenvalues for the first 
three principal components were such as to explain generally over ninety percent 
of the variation in the data. The eigenvectors have also been calculated for 
each image and for the fourteen band composite obtained by combining the two 
•even band images: in this case the first seven bands consisted of those from 

one image and the latter seven bands (in the same sequence) contained digital 


radiance data from the aacond image. Tha avaluatad aitan vector* are shown in 
Table 2. 

FI*. 4 show* a plot of tha eigenvectors for tha firat principal component 
for each data, for axle rotationa par formed to oaphaalsa land uae type* in each 
training area and for tha whole imatre extract. Tha aifanvaetora for tha aacond 
and third principal component* ahow leaa diffaranca between land uaa type* 
categorized by tha trainin« aitaa for tha combined or for tha single-date 
images. However, in each caaa, there ia a aubatantial difference between each 
of tha eigenvector* for a fiven principal component and the eigenvector for the 
same principal component for the whole image. Thia demonatrate* the eepar- 
ability of aach land ua* type using this technique. 

Fip. 5 shows the percenter* of variation in the data explained by each 
principal component for the fourteen band image. This is the only image so 
illustrated, for reasons of space. Fip. 6 shows the loadings (scaled to the 
rang* 1.0) of each of the components (feature vectors) in the eigenvector for 
each principal component calculated to emphasize the four features contained in 
training areas considered for the two images (fourteen band) composite. The 
forest area clearly has considerable shadow and a low albedo and thus has a 
lower eigenvector (first principal component) than that of the dense urban area, 
which is less than that of the airport, which in turn is less than that of the 
downtown area. Note that the loadings of the components (measurement vectors) 
for each of the training sites shows a progressive change from forest to down- 
town area. It is suggested that both shadow and the vegetative component con- 
tained in the training area contribute to the differences in loadings. It is 
noted that the loadings for the bandpasses of the March 1964 images generally 
exceed (in modulus) those for the November 1962 image. The seasonal vegeta- 
tional change and the change in shadow component due to the change in solar 
azimuth and zenith angles at the time of satellite overpass will also be 



important In contributing to not* information in tht coeond iM|t which was not 
present In tha flrat. This explains why the first principal component calcu- 
lated for tha composite fourteen band image was superior in separating the 
features contained within the training areas identified earlier on and shown in 
Fig. 2. 

It is suggested that either a canonical variate formed from the evaluated 
eigenvectors obtained by training on a class containing a land use category or 
possibly the principal components of ratioed images might be worth investigating 
as bstter land use type discriminators. Clearly, these possibilities need to be 
investigated and the Intercalibration between the digital radiance data from the 
Landsat 4 and 5 thematic sensors needs to be taken into cccount (e.g. Duggin 
1965, Likens and Urlgley 1965. Palmer 1984). 

CPNCLVSIPy 

A feasibility study Involving the principal components analysis performed 
by training upon selected urban land use types to make them most separable from 
other imaged features has been described for single-date and for registered 
multidate images. It appears that separability is enhanced by using the 
increased spectral radiance information contained in superimpesed images 
obtained on different dates. It appears probable that shadow and vegetation 
vigor, as well as changes in these parameters are importoat in controlling 
feature separability. It appears that TH bands 4 and 6 contain information which 
best discriminates forest from other areas and which stay therefore be of use in 
urban forest studies. Further work on a larger suite of images, containing a 


larger variety of urban forest land use patterns is underway. 
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TABLE CAPTIOH 


Tabic 1. Mean radiance values for each training site (land use category) for 
each inage studied. 

Table 2. Calculated eigenvectors for the principal components obtained by 
training on selected features, to so rotate the axes as to enhance the 
separability of those features selected and by calculating the principal 
components axis rotation for the whole image. Values are shown for the first 
three principal components for each of the two images and for the fourteen band 


superimposed image pair 


FIGURE CAPTIONS 


Fig. 1 

Fig. 2 
Fig. 3 

Fig. 

Fig. 5 

Fig. 6 


Single band extracts of images euper imposed Lendsat thematic mapper 
(512 pixels x 512 lines) obtained on two dates; 2 November 1982 
(Landsat 4) and 24 March 1984 (Landsat 5). The images are contrast 

enhanced by using a linear stretch. 

Training areas used in 512 x 512 pixel image. 

Principal components images obtained for the whole 512 pixel x 512 
line extract of the fourteen band composite of the superimposed 
images. The training area (extracts) of the few selected cover types 
are also shown. 

. Plot of the eigenvectors as a function of training site (Fig. 2) and 
for the whole 512 x 512 pixel image extraction for the first principal 
component for each of the images and for the composite of the two 
superimposed images. 

Plot of the percentage of variation in the data explained by each of 
the principal components for the fourteen band (composite) image. In 
each case the axes have been so rotated to as facilitate the 
extraction of the selected feature. 

Loadings (scaled to fall between 1.0) of the feature vectors in 
calculating the eigenvectors for the first three principal components 
when the axes are so rotated as to emphasize each of the selected 


features. 
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Table 1 


TW 


band 

Foraat 

Airport 

Downtown 

Danse 

Urban 


Nov. 62 

March 64 

Nov. 82 March 84 

Nov. 82 

March 84 

Nov. 82 March 84 

1 

58.7 

86.3 

72.1 

99.7 

70.9 

99.9 

67.8 

95.3 

2 

23.6 

33.6 

30.5 

42.0 

29.4 

42.0 

27.7 

38.7 

3 

24.6 

36.0 

31.4 

46.4 

30.3 

45.9 

26.5 

42.5 

4 

37.3 

49.5 

35.2 

47.6 

33.6 

50.0 

31.6 

42.8 

5 

40.8 

76.6 

48.0 

75.9 

42.9 

74.4 

41.9 

70.6 

6 

113.5 

112.1 

116.0 

121.9 

116.5 

120.0 

117.2 

124.0 

7 

16.0 

33.2 

22.8 

36.2 

21.3 

37.3 

21.6 

37.7 

6/4 

3.04 

2.26 

3.29 

2.56 

3.46 

2.40 

3.71 

2.90 
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Table 2. 

PC 

DATE 

FOREST 

DENSE 

AIRPORT 

DOWNTOWN 

WHOLE 




URBAN 



IMAGE 

1 

Nov 82 

-48.11 

-35.94 

-21.64 

65.1 

85.52 


March 64 

+42.25 

51.85 

74.03 

132.24 

126.77 


Nov 82 + 
March 84 

-50.. 03 

-28.41 

+24.75 

+142.56 

152.84 

2 

Nov 82 

130.50 

137.02 

150.79 

102.78 

50.04 


March 84 

169.29 

168.82 

172.94 

77.01 

89.79 


Nov. 82 + 
March 84 

+200.95 

206.11 

212.01 

110.84 

108.06 

3 

Nov 82 

23.78 

-0.34 

10.85 

66.94 

26.25 


March 64 

-18.24 

-4.80 

-22.75 

44.67 

13.99 


Nov 82 + 
March 84 

55.22 

74.21 

57.40 

81.31 

23.40 
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